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Abstract 



We perform a detailed study of the decays of the lighter top and bottom squarks 
(ti and 61) in the Minimal Supersymmetric Standard Model (MSSM). We show 
that the decays into Higgs or gauge bosons, i.e. t% — ► b\ + {H + or W + ) and 
bi — > t\ + (H~ or W~), can be dominant in a wide range of the MSSM parameters 
due to the large Yukawa couplings and mixings of t and b. Compared to the decays 
into fermions, such as t\ — ► t + (xi or g) and t\ — > b + xj, these bosonic decay modes 
can have significantly different decay distributions. We also show that the effect of 
the supersymmetric QCD running of the quark and squark parameters on the t\ 
and 61 decay branching ratios is quite dramatic. These could have an important 
impact on the search for t\ and b\ and the determination of the MSSM parameters 
at future colliders. 



In the Minimal Supersymmetric Standard Model (MSSM) [I] supersymmetric (SUSY) 
partners of all the Standard Model (SM) particles are introduced. In order to solve the 
hierarchy, fine-tuning and naturalness problems, the SUSY particle masses have to be less 
than 0(1 TeV). Hence the discovery of all SUSY partners and the study of their properties 
are essential for testing the MSSM. Future colliders, such as the Large Hadron Collider 
(LHC), the upgraded Tevatron, e + e~ linear colliders, and colliders will push the 

discovery reach for the SUSY particles up to the TeV mass range and allow for precise 
measurement of the MSSM parameters. 

In this article we perform a phenomenological study concerning the search for the SUSY 
partners of the top (t) and bottom (b) quarks (i.e. stops (t) and sbottoms (6)). These 
particles may have properties which are very different from those of the squarks of the 
other two generations due to the large top and bottom Yukawa couplings. Production and 
decays of t and b were studied in @, || f|, ||. Stops (sbottoms) have two mass eigenstates 
ti,2 (^1,2) with mi x < mi 2 (mg < wg 2 ). Here we focus on tx and b\. Like other squarks, 
they can decay into fermions, i. e. a quark plus a gluino (g), neutralino (xl) or chargino 
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h -> t + (g or Xi ) , h -> b + (g or xii , , , , 

h -> b + x} > 61 -> * + Xj , 

with i = I, ...,4 and j = 1,2. In addition, they can also decay into bosons 0, |[], i. e. a 
squark plus a Higgs or gauge boson: 

i x -> 61 + (#+ or W+) , 61 -> ti + (F- or W) . (2) 

The decays of Eq. (Q) are possible in case the mass difference between tx and b\ is suffi- 
ciently large. 

In the present article we extend the analysis of 0, [|. We point out that the bosonic 
decays of tx and 61 of Eq. (0) can be dominant in a large region of the MSSM parameter 
space due to large top and bottom Yukawa couplings and large t and b mixing parameters. 
This dominance of the bosonic decays over the conventional fermionic decays of Eq. ([I]) 
could have an important impact on searches for tx and 61 at future colliders. An analogous 
study for t 2 , b 2 , f 2 , and v T was performed in || |[, where f 2 and v T are the sleptons of the 
third generation. In addition we show that the effect of SUSY-QCD running of the quark 
and squark parameters @ on the decay branching ratios of tx and bx is quite dramatic. 

First we summarize the MSSM parameters in our analysis. In the MSSM the squark 
sector is specified by the mass matrix in the basis (q L , q R ) with q — t or b [R Q 



with 



a q m q m~ R 



m \ = M 2 Q + m 2 z cos 2(3 (II L - e q sin 2 6 w )+m 2 qJ (4) 
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m 2 g R = M 2 {(j t)] + m% cos 2(3 e q sin 2 9 W + m 2 , (5) 

{ {A t - n cot 13) mt (q = t) 
q q " \ {A b -fitan/3)m b (q = b). { ) 

Here If is the third component of the weak isospin and e q the electric charge of the quark 
q. Mq q f) and A t)b are soft SUSY-breaking parameters, \x is the higgsino mass parameter, 
and tan/3 = v 2 jv\ with v x {v 2 ) being the vacuum expectation value of the Higgs field H® 
(H$). We treat Mq ^g and A tjb as free parameters since the ratios Mq/Mq, M^/Mq and 
A t /A b are highly mo del- dependent. By diagonalizing the matrix (|3|) one gets the mass 
eigenstates q x = q L cos 9g + q R sin 9g, q 2 = —q^ sin 9g + q R cos 9g with the mass eigenvalues 
m gi5 m Q2 ( m 9i < m ft>) an d the mixing angle 9g (— | < 9g < |). As can be seen, sizable 
mixing effects can be expected in the stop sector due to the large top quark mass. Likewise, 
bL-bR mixing may be important for large tan/3. 

The properties of the charginos xt (* — 1> 2; < m ^±) an d neutralinos x° (k = 1, 4; 
m^o < ... < m^o) are determined by the parameters M, M', // and tan /?, where M and M' 
are the SU(2) and U(l) gaugino mass parameters, respectively. Assuming gaugino mass 
unification we take M' = (5/3) tan 2 OwM and m g = (a s (rrig) / 'a 2 )M with m g being the 
gluino mass. The masses and couplings of the Higgs bosons h°, H°, A and H ± , including 
leading radiative corrections, are fixed by m^, tan/3, //, m t , m b , Mq, Mfj, Mg, A t , and 
A b . H° (h°) and A are the heavier (lighter) CP-even and CP-odd neutral Higgs bosons, 
respectively. For the radiative corrections to the h° and H° masses and their mixing angle 
a we use the formulae of Ref . |10| ; for those to m#+ we follow Ref . JTTJ] []. 



The widths of the squark decays into Higgs and gauge bosons {H and W ) are given 
by 1: 



rfo-fcir*) = 7 ^g 2 + , (7) 



167T m| x 
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lo7r m^m^ 

Here (qi,q[) = (ti, &i) or = ^(m^, m|/, m 2 ^) is the usual kinematic factor, 

K,(x,y, z) = (x 2 + y 2 + z 2 — 2xy — 2xz — 2yz) 1 / 2 . Notice an extra factor k 2 /m^ for the 
gauge boson mode. The G H + denotes the squark coupling to H + and C w + that to W + . 
Their complete expressions, as well as the widths of the fermionic modes, are given in 

EE- 

From Eqs.(3-6) we see that mi 12 ~ Mq ^ and m^ l2 ~ Mq^ i n case ^Q,u,b are l ar g e 
relatively to the other parameters. In this case, for Mfj > Mq 3> Mg (Mg > Mq 3> M^) 
we have ^> m~ bi (m^ 3> m.^ ) , which may allow the bosonic decays of Eq. (Q) . Hence 
in this article we consider two different patterns of the squark mass spectrum: mi x ^> m bl 
with (fi,6i) ~ (t L ,b R ) for M > Mq > Mg, and > m fl with (f^fti) ~ (*«,&£,) for 
1 Notice that fil], 111) have a sign convention for the parameter fj, opposite to the one used here. 
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Mjj 3> Mq ^> Mjj. Note here that in the former (latter) pattern the condition My 3> Mq 

(Mfj 3> Mq) ensures t\ ~ t L (b\ ~ which eventually enhances the bosonic decays of 
Eq.(^) as we will see below. Thus the bosonic decays considered here are basically the 
decays of into 6r and b \ into {r. 

The leading terms of the squark couplings to if ± and W ± are given by 

Gh+ = G(tibiH ± ) ~ h t (fM sin/3 + A t cos/3) sin%cos 0g 

+ /i&(//cos/3 + A 6 sin/3) cos^sin^, (9) 

C w + = C^tiSiW^) ~ -^cos%cos0 g . (10) 

V 2 

The Yukawa couplings /i tj b are given as 

h t = gm t /(V2m w sin/5), h b = gm b /(V2m w cos (3). (11) 

The Higgs bosons if ± couple mainly to Qlq'r combinations. These couplings are propor- 
tional to the Yukawa couplings ht,b and the squark mixing parameters A tj b and /x, as can 
be seen in Eq.(|9|). Hence the widths of the squark decays into H ± may be large for large 
A tt b and /i. Note here that the squark mixing angles 9g themselves depend on A t b , /i and 
tan j3. In contrast, the gauge bosons W ± couple only to <7l<?l, which results in suppression 
of the decays into W . However, this suppression is largely compensated by the extra 
factor k 2 /m^, in Eq.(|8|) which is very large for — m% ^> mw- (This factor stems from 
the contribution of the longitudinally polarized gauge boson radiation (qi — > q[W^).) 
Hence the widths of the squark decays into W ± may be large for a sizable q' L — q' R mixing 
term a q /m q i. On the other hand, the fermionic decays are not enhanced for large A t ^ 
and /i. Therefore the branching ratios of the bosonic decays of Eq.(Q) are expected to be 
large for large A t b and /i if the gluino mode is kinematically forbidden. As for the tan f3 
dependence, we expect that the branching ratio B(ti — > b\ + (H + ,W + )) increases with 
increasing tan/3 while B(b\ — > t\ + (H~ , W~)) is rather insensitive to tan/3. The reason 
for this is as follows: 

(i)fi decay; 

As (ii,&i) ~ (^l,&r) in this decay, the coupling Gh+ = G{tib\H + ) ~ G(ik&Rii + ) oc h b oc 
tan/?. Hence we expect that the width T(t\ — > b\H + ) ( and hence B{t\ — > b\H + )) in- 
creases with increasing tan/?. The coupling Cw+ = C{tibiW + ) increases with increasing 
bi — bn mixing term a^mfe which increases with increasing tan/3. Hence we expect an 
increase of B{ti — > b\W + ) as tan/3 increases, 
decay; 

As (£i,&i) ~ (tR,bi) in this decay, the coupling Gh- = G{t\biH~) ~ G{tRbLH~) oc h t . 
Hence we expect that B(b\ — > t\H~) is rather insensitive to tan/3 ( for tan/3 > 3). The 
coupling Cw- = C(tibiW~) increases with increasing — t R mixing term a t m t which 
is fairly insensitive to tan/3 (for tan/3 > 3 ). Hence we expect a mild dependence of 
B(b! tiW~) on tan/3. 
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We now turn to the numerical analysis of the t\ and 61 decay branching ratios. We 
calculate the widths of all possibly important two-body decay modes of Eqs. (|IJ) and (0). 
Three-body decays are negligible in this study. The widths of the t\ and b\ decays into 
receive large SUSY-QCD corrections for large tan/3 in the on-shell renormalization 



scheme |12| : the 0(a s ) correction terms are often comparable to or even larger than the 
lowest order ones. Such large corrections make the perturbation calculation of the decay 
widths unreliable. In general this problem shows up in calculating rates of processes 
involving the bottom- Yukawa-coupling {h b ) for large tan/3 in the MSSM [|7]]. In Ref. |7j 
it is pointed out that this problem can be solved by carefully defining the relevant tree- 
level couplings in terms of appropriate running parameters and on-shell squark mixing 
angles 9 q . Following Ref. f7j, we calculate the tree-level widths of the ti and bi decays, 
Eqs.([TJ) and (|2|), by using the corresponding tree-level couplings defined in terms of the 
SUSY-QCD running parameters m q (Q) and A q (Q) (with the renormalization scale Q 
taken as the on-shell (pole) mass of the decaying squark), and the on-shell squark mixing 
angles 9 q . For the kinematics, i.e. for the phase space factor n/m^ (for q~i decay) the 
on-shell masses are used. We call the widths thus obtained as 'renormalization group 
(RG) improved tree- level widths'. Here we adopt the notations and conventions of J7J. 
Our input parameters are all on-shell ones except A b which is a running one, i.e. they are 
M t , M b , Mq(t), Mjy, M^, A t , A b (Q), /i, tan/3, m^, and M, with the renormaliztion scale Q 
at the on-shell (OS) mass of the decaying squark m^os- M tjb are the on-shell (pole) masses 
of the t,b quarks. Mg(q) is the on-shell Mq for the q sector. Note that Mg(t) is different 
from Mq(6) by finite QCD corrections [|7|. The procedure for obtaining all necessary 
on-shell and _Di?-running parameters of quarks and squarks (i.e. on-shell (pole) m qi2 , 
on-shell 9 q , and SUSY-QCD running (m q , A t , Mq, Mjj, M^,, m qi 2 )) is described in detail 
in @. We use this procedure. For the SM parameters we take M(=175GeV, Mf>=5GeV, 
mz=91.2GeV, sin 2 9w = 0.23, mw = mz cos 9 w, a (rriz) = 1/129, and a s {mz) = 0.12 
[with the one-loop running a s {Q) = 127r/((33 — 2nj) ln(Q 2 / A^ f )) , Uf being the number 
of quark flavors; only for the calculation of the SM running quark mass m q (Q)sM from 
the two- loop renormalization group equations we use the two-loop running a s (Q) as in 
0]. In order not to vary too many parameters we choose Mq{t) = \Mfy = |M^, (Mq{1) 

= \Mfj = \Mf ) ) for t\ (61) decays, and A t = A b (Q) = A (where Q=m qi os f° r <li decay) 
for simplicity. Moreover, we fix M=400GeV (i.e. mg=1065GeV) and mA=150GeV. Thus 
we have Mq(£), A, /1 and tan/3 as free parameters. 

In the plots we impose the following conditions: 



(i) m-± > 100 GeV, m h o > 110 GeV, m fi i iOS > m^o, m^o > 90 GeV, 

(ii) 5p (t—b) < 0.0012 ]13| using the formula of |14| (the constraint from electroweak 5p 
bound on t and 6), and 

(iii) A*(Q) < 3 (M^+M^+m^), and Al(Q) < 3 (M%(Q) +M|(Q) +m 2 Hl ) with 
m H 1 — ( m A + m z) s i n2 P ~ \ m i> m H 2 = ( m A + m l) cos2 P ~ I m l; an d Q ~ Mq 
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(the approximate necessary condition for the tree-level vacuum stability [p~5j| ) . 

Conditions (i), along with m- g = 1065 GeV, ensure that we satisfy the experimental 
bounds on xf, Xi, h°, ti, bi and g from LEP2 |16| and Tevatron |[7|]. Note that m^o > 90 
GeV is imposed in order to evade the experimental bounds on g . Conditions (ii) and 
(hi) constrain the t and h mixings significantly. We note that the experimental data for 
the b — > S7 decay give rather strong constraints |Tj| on the SUSY and Higgs parameters 



within the minimal supergravity model, especially for large tan (3. However, we do not 
impose this constraint since it strongly depends on the detailed properties of the squarks, 
including the generation-mixing. 

In Fig.l we plot in the A-fi plane the contours of the t\ decay branching ratios of the 
Higgs boson mode B(t\ — > b\ + H + ), the gauge boson mode B(t\ — > b\ + W + ), and the 
total bosonic modes B(i x -> h + (H+, W+)) = Bfc -> h + F+) + Bfo -> 6 X + 
at the RG-improved tree-level for tan/3=30 and Mg(t)=600GeV. We show also those of 

5(tx — ► b\ + (H + , W + )) at the naive (unimproved) tree-level, where all input parameters 
are bare ones, i.e. they are m t (=175GeV), m 6 (=5GeV), Mq = \M = \M b (=600GeV), 
A t = A b = A, ft, tan/3 (=30), m A (=150GeV), and M(=400GeV) (see Eqs.(|)-(D). We see 
that the t± decays into bosons are dominant in a large region of the A-fi plane, especially 
for large \A\ and/or as we expected. Comparing Fig.l.c with Fig.l.d we find that the 
effect of running of the quark and squark parameters (m q (Q), A g (Q), Mq ^ f)(Q)) is quite 
dramatic. Note here that the symmetry under A — > —A and/or /j, — > — //, which holds at 
the naive tree-level, is strongly broken at the RG-improved tree-level. This is mainly due 
to the running effect steming from the gluino loops 0. 

In Fig. 2 we show the contours of the b x decay branching ratios B(b\ — > ti+H~),B(bi — > 
t\ + W~), and B(b\ — > t\ + (H~,W~)) at the RG-improved tree-level for tan/3=8 and 
30, and Mg(t)=600GeV. The results are very similar to those for the t\ decays. We have 
found again a dramatic effect of the running of the parameters. 

In Fig. 3 we show the individual branching ratios of the t x and b\ decays as a function of 
tan/3 for (A, /i, Mg(t)) = (-800, -700, 600)GeV and (800, 800, 600)GeV, respectively. One 
can see that the branching ratios of the t\ decays into bosons increase with increasing 
tan/3 and become dominant for large tan/3 ( > 20), while the b\ decays into bosons are 
dominant in the entire range of tan/3 shown, as expected. As already explained, these 
tan (3 dependences of the i\ and b\ decays come from the increase of (h b , ctbrrib) with tan f3 
and the mild dependence of {h t ,a t m t ) on tan/3, respectively. 

In Fig. 4 we show the Mg(t) dependence of the t\ and b\ decay branching ratios for 
(A(GeV), / u(GeV),tan/3) = (-800, -700, 30) and (800, 800, 30), respectively. In these cases 
we have (m^o, m-+, m^) = (198, 393, 1065)GeVand (198, 394, 1065)GeV, respectively. One 

can see that the bosonic modes dominate the t\ and b\ decays in a wide range of Mg(t). 
For the &i decays we have obtained a similar result for (A(GeV), //(GeV), tan/3) = (800, 
-800, 8). (Note that the decay into a gluino becomes dominant above its threshold.) 
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One can also see that the branching ratios of the bosonic decays decrease with increasing 
Mqit) which is roughly equal to the mass of the decaying squark q x . This comes from the 
fact that in the large Mg(i)(~ m^) limit the decay widths of the bosonic and fermionic 
modes are proportional to m^ 1 and , respectively. 

We find that the dominance of the bosonic modes is fairly insensitive to the choice 
of the values of ttia, M, and the ratio Af ) (Q)/A t . The decays into H ± are kinemati- 
cally suppressed for large m^. However, the remaining gauge boson mode can still be 
dominant. For example, for Mg(t)=600GeV, A=/j=800GeV, and tan/?=30, we have 

(m t ~ iOS , m^ lOS , rn^o, m-+) = (347, 602, 198, 394)GeV with S(6i -> fi#-)=(48, 37, 13, 0)% 

and B{b x -> h W~)=(43, 51, 71, 82)% for m A =(120, 200, 240, >250)GeV. We have also 
checked that our results do not change significantly if we take smaller values of M, as long 
as decays into a gluino are kinematically forbidden. As for the sensitivity of our results 
to the ratios Mjj f}/Mn(t), we have found that for small tan/3 and large \A\ and |//| the 
bosonic modes can dominate the b± decay even in case Mg(t) = My = M^; e.g. for 
tan/3 = 8 and Mq(t) = Mjj = Mfj = 600GeV, we have obtained results similar to Fig.2.b 

and Fig.2.c. For example we have B(bt -> t x +W~) = 60%, S(6i -> ii+(H~, W~)) = 90% 
and (mi lOS ,mi 20S ,mj )iOS ,mj )20S ,m H +) = (388, 792, 570, 638, 166)GeV for A=1200GeV, 
/i=-1300GeV, tan/3 = 8, and M Q (t) = M = M b = 600GeV. For these parameters the 
tL-tn mixing term is large and the &l-6r mixing term is small. This leads to a large 
mass-splitting of tyti and a small one of &1-&2, which results in sufficient phase spaces for 
the bosonic decays of b\. 

In a complete analysis one would have to calculate the full SUSY-QCD one-loop cor- 
rections to the widths of the ti and bi decays. However, we expect that these corrections 
will not invalidate the dominance of the bosonic modes in the ti and b\ decays in a sig- 
nificant portion of the MSSM parameter space as shown in our analysis. The calculation 
of the full corrections is beyond the scope of the present paper. 

Now we discuss the signatures of the t\ and b\ decays. We compare the signals of the 
decays into bosons (Eq. (|)) with those of the decays into fermions (Eq. (P). In principle, 
the final states of the bosonic decays can also be generated from fermionic decays. For 
example, the final particles of the decay chain 

k + (H + or W + ) -> (&x?) + (qf) (12) 

are the same as those of 

h - b + xt,2 ~+ b + ((H + or W+) + *?) - b + (qq'xi). (13) 



Nevertheless, the decay distributions of the two processes (|i2|) and (0) are in general 
different from each other due to the different intermediate states. For example, the b 



in the chain ([12]) tends to be softer than the b in (0). A similar argument holds for 
the quark pairs qq' in the decay chains. Moreover, the distribution of the missing energy- 
momentum carried by Xi could be significantly different in (12) and (|13D since it is emitted 
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from a different sparticle. Hence the possible dominance of the bosonic decays over the 
conventional fermionic decays could have an important impact on the search for the ti 
and b\, and on the measurement of the MSSM parameters. Therefore, the effects of the 
bosonic decays should be included in the Monte Carlo studies of the t\ and bi decays. 

In conclusion, we have shown that the t\ and b\ decays into Higgs or gauge bosons, 
such as t\ — *■ b% + (H + or W + ), can be dominant in a fairly wide MSSM parameter region 
with large mass difference between ti and b\, large \A ttb \ and/or |//|, and large rrig (and 
large tan/3 for the t\ decay) due to the large Yukawa couplings and mixings of t and 
b. Compared to the fermionic decays, such as t\ — > b + xt,2i these bosonic decays can 
have significantly different decay distributions. We have also shown that the effect of the 
SUSY-QCD running of the quark and squark parameters on the t\ and b\ decays is quite 
dramatic. These could have an important impact on the searches for t\ and &i and on the 
determination of the MSSM parameters at future colliders. 
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Figure Captions 



Figure 1: Contours of branching ratios of t\ decays at the RG-improved tree-level in the 
A-/x plane for tan/3 = 30, and Mg(t) = \M d = \M & = 600GeV; (a)B(t 1 -> bi + H + ), 

(b)S(fi -»• 61 + W + ), and (c)S(ti -> 61 + (ff + , W + )). The regions outside of the dashed 
loops are excluded by the kinematics and/or the conditions (i) to (iii) given in the text. 
Contours of the corresponding branching ratio B(ti — > 61 + (ff + , W + )) at the naive tree- 
level are shown in Fig.d. 

Figure 2: Contours of branching ratios of b\ decays at the RG-improved tree-level in 
the A- fx plane for tan/3= 8 (a-c) and 30 (d-f), and Mg(t) = §A% = fAfg = 600GeV; 

(a,d)S(6i -> ti + if"), (b,e)S(6i -»• 1 1 + W~), and (c,f)S(6i -> 1 1 + (ff-,Vy-)). The 
regions outside of the dashed loops are excluded by the kinematics and/or the conditions 
(i) to (iii) given in the text. 

Figure 3: tan (3 dependence of t\ (a) and b\ (b) decay branching ratios for (A, fx, Mq(1)) 
= (-800, -700, 600)GeV and (800, 800, 600)GeV, respectively. "61 + H + /W + " and % + 
H~/W~" refer to the sum of the Higgs and gauge boson modes. The gray areas are 
excluded by the conditions (i) to (iii) given in the text. 

Figure 4: Mg(t) dependence of t\ (a) and b\ (b) decay branching ratios for (A(GeV), 

At(GeV), tan/3) = (-800, -700, 30) and (800, 800, 30), respectively. "61 + H + /W + " and 
"ti + H /W~" refer to the sum of the Higgs and gauge boson modes. The gray areas are 
excluded by the conditions (i) to (iii) given in the text. 
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Fig.l 
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Fig.4 
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